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Abstract: In the Friedmann cosmology the deceleration of the expansion q plays a fun- 
damental role. We derive the deceleration as a function of redshift q{z) in two scenarios: 
ACDM model and modified Chaplygin gas {MCG) model. The function for the MCG 
model is then fitted to the cosmological data in order to obtain the cosmological param- 
eters that minimize x 2 - We use the Fisher matrix to construct the covariance matrix of 
our parameters and reconstruct the q(z) function. We use Supernovae la, WMAP5 and 
BAO measurements to obtain the observational constraints. We determined the present 
acceleration as qo = —0.60 ± 0.12 for the MCG model using the Constitution dataset of 
SNela and BAO, and qo = —0.63 ± 0.17 for the Union dataset and BAO. The transition 
redshift from deceleration to acceleration was found to be around 0.6 for both datasets. We 
have also determined the dark energy parameter for the MCG model: Oxo = 0.834 ±0.028 
for the Constitution dataset and Qxo = 0.854 ± 0.036 using the Union dataset. 
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1. Introduction 

During the last decade, the observation of type la supernovae (SNela) and the cosmic mi- 
crowave background radiation (CMBR) permitted the determination of the cosmological 
parameters with ever increasing precision. The reported results of the seven year analysis 
of WMAP |] are Q A = 0.734 ± 0.029 and 9. m = 0.266 ± 0.029; these values were obtained 
assuming a flat geometry {VL^o = 0). From measured luminosity distances to SNela, Riess 
et al. 0, |3| determined the redshift of the transition from decelerated to accerelerated ex- 
pansion to be zt = 0.46 ±0.13; this value was obtained assuming a linear expansion for the 
deceleration parameter, that is, q(z) = qo + q\z. The problem with this linear expansion 
is that it works well for small redshifts but the transition redshift is not so small. Using 
another parametrization for q(z) Shapiro and Turner Q concluded that the present SNela 
data cannot rule out the possibility that the universe has been decelerating since z = 0.3. 
In other references ( || , || ) it was shown that the value of the transition redshift depends 
on the adopted parametrization for q(z), as well as on the data sample. In particular, a 
parametrization that has the redshift transition as free parameter has been presented in 
reference [^j . This parametrization could be used to study the kinematics of the expansion 
regardless of the matter content of the Universe. However, the statistical properties of 
data are still not good enough to produce strong constraints. In general, the parametriza- 
tions are helpful by their phenomenological properties, since they can serve to study the 
accelerated expansion in different contexts, such as the structure formation, that would be 
difficult to study into a fundamental theory. 
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Among all possible candidates to explain the accelerated expansion, the Chaplygin gas 
is a strong candidate; it is the best known proposal of a unification of dark matter {DM) 
with dark energy (DE) into a single fluid. The ideia is: an equation of state (EoS) leads 
to a component which behaves as dust at early stage and as cosmological constant at later 
stage. Following this idea it was considered the so called generalized Chaplygin gas This 
model has been analyzed many times in the literature; see for example ||. Subsequently 
this model has been modified to include an initial phase of radiation, and is called the 
modified Chaplygin gas (MCG). From the theoretical point of view, this scenario can 
also be restated as a Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmological model 
containing a scalar field <j) with its self-interacting potential [llj. 

In this paper we have as main aim to study the MCG model using several data to con- 
straint as much as possible the parameters of the model. In the literature the generalized 
Chaplygin gas has been studied in the context of statefinder diagnostic, stochastic gravita- 
tional waves, observational constraints using gamma ray bursts, strong lensing, Supernovae 
la, etc.; see [11 1 for references. We studied the deceleration parameter in the MCG and 



used recent observational data to constraint its free parameters. We used two sets of data 
of type la Supernovae: the Constitution set and the Union set. We also use the CMB and 
BAO data. In order to make comparisons, we derive the function q(z) for the standard 
cosmological model, that is the ACDM model. The derivation was performed with the 
minimal assumptions: GR theory is valid and the universe is homogeneous and isotropic, 
that is, the cosmological metric is the FLRW one. No assumption was made about the 
spatial curvature of the universe, however the analysis on the observational limits of the 
free parameters was restricted to flat models. We leave for a future paper the models with 
curvature. 

The outline of this paper is as follows. In section 2 we present the equations of the 
models, ACDM and MCG, and determine the Hubble parameter and the deceleration 
parameter respectively. In section 3 we describe the method used to obtain the confidence 
regions and the reconstruction of the function q(z). Finally, section 4 is devoted to the 
discussion of our results. 

2. Equations of our Models 
2.1 ACDM Model 

Let's assume that the universe is described by the FRW metric and the energy content is a 
pressureless fluid and a cosmological constant. The first Friedmann equation for the scale 
factor a(t) reads 



^ 2 = _ 



a 2 _ SirGp A k 
72 ~^S~ + 3~~2 



which can be written as 



with the definitions 



n m + n A + n k = 1 (2.2) 
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The Bianchi identity for a pressureless fluid can be integrated to give 



pa 3 = p al 



(2.3) 



where the subscript zero denotes present values of the quantities. Inserting eq.( p.3| ) into 
eq. (12. II) and writing — k as 



k = H 2 a 2 Q k = H^a z o n ko 



2 n 2, 



we obtain 



H a 



(2.4) 



In terms of the redshift z = ^ — 1, eq.(2.4) reads 



(1 + *) 
where 



,2_ H 2 a 2 E 2 {z)=H 2 a 2 iz)E 2 {z) (25) 



^ 2 (z) = n m0 (i + z) 3 + n k0 (i + z) 2 + n A0 (2.6) 



Using eq.(2.5) we can relate % with d = 



da 

dz 



«o 



a(z) 



da 



(l + z) 2 l + z H (l + z)E(z) dt 
from which we infer the relation between dt and dz: 



(2.7) 



dz 
~dl 



-H (l + z)E(z) 



Now, using eq.(|2.7D and calculating d, we obtain for the deceleration parameter 



H 2 c 



the expression 
q(z) = 



—-(l+z) -^A 



[fi ro0 (l + *) 3 + «*o(l + *) 2 + «Ao] 1 (2. 



For z » 1 we see that q(z) = 0.5 as expected. The transition redshift, where q(zt) = 0, 
is given by 
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for any value of the curvature k. 

In Figure 1 we reconstruct the evolution of q(z) using the Constitution dataset and 
Union dataset. From this figure we find that both datasets predict similar transition 
redshifts, z t ~ 0.5. Our reconstruction was made with a la confidence level. 



l.Or 



l.Or 
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Figure 1: We show the observational constraints for the AC DM model using only SNela with 
fibo = 0.042 ± 0.027. Left panel: Constitution set; right panel: Union set. 



2.2 Modified Chaplygin gas (MCG) model without cosmological constant 

The generalized Chaplygin gas model has been proposed as a source term in Einstein's field 
equations in order to unify the concepts of cold dark matter and dark energy [b|, [llj . 
We will consider the energy content of the universe as a fluid that behaves like a perfect 
fluid of non-zero pressure at early times, like a pressureless fluid at intermediate times and 
like dark energy at present. 

The equation of state of the MCG is given by 

p= Bp-^, (2.10) 

where A,B and a are non-negative constants. For B = we have the pure generalized 
Chaplygin gas and for A = a perfect fluid. The MCG behaves as radiation (when 
B = 1 /3) or dust-like matter (when B = 0) at early stage, while as a cosmological constant 
at later stage. On the other hand, the Bianchi identity 

3-(p + p) + p = (2.11) 
a 

yields after integration (see |l5|| for details): 
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Pifl) = Po 

where 



n x + (i-n x ) (^) 3R 



l/(a+l) 

(2.12) 



R=(B + l)(a + l) 
Po is the present energy density, and we define the dimensionless parameter 

n x = tt • (2-13) 

The same analysis of the previous section can be done, leading to the following expres- 
sion for the deceleration parameter: 

= ^ [n* + (i - n x )(i + z f^]- al(l+a) 



{(i - n x )(w + i)(i + zf R - 2n x } (2.14) 



where 



E 2 (z) = [n x + (1 - Q x )(l + z) 3 «] 1 /(i+«) + Ofc0 (l + z) 2 (2.15) 
The transition redshift for the MCG model is given by: 

< 1 + *'> 3 *- (33 + lKl-» ri - (2 ' 16) 

The dimensionless parameter represents the fraction of dark energy in the content 
of the universe, thus taking the value £l x ~ 0.7. If we chose B = 1/3 in order to describe 
the evolution of the universe since the radiation era, and consider < a < 0.5, which is 
required from thermodynamical considerations (see Jl6|]), we find z% ~ 0.2. If we consider a 
pure Chaplygin gas (B = 0) and ft x « 0.7, the value Zt = 0.46 can be achieved if a « 0.4. 



This result is compatible with the analysis done in [17] for a Chaplygin gas in the flat 
(k = 0) case. 

Now, in order to carry out our analysis of observational constraints, we consider as 
components of the Universe: baryons plus MCG. Thus, in the flat case, the Hubble 
parameter is given by: 



e 2 (z) = f^M) 2 = n w (i + z f + (i - n w )[(i - n xo )(i + z) m + Oxo)] 1/(1+Q) (2.17) 

This expresion will be used in our analysis of observational constraints in the following 
section. As can be seen, we have four free parameters (fi&Qj ^xo, ot, B). 
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3. Some Observational Constraints 



In the present section we consider some observational constraints of SNela, BAO and CMB 
for our models. In this context it is important to consider the comoving distance to an 
object at redshift z, 

where we consider only the fiat case; E{z) = H(z)/Hq is given for the ACDM model by 
equation (2.6) and for the MCG model by equation (2.17). Using the equation above, the 
luminosity distance in the flat case is then given by di, = (1 + z)r(z). 

3.1 Constraints from Supernovae Data 

The supernovae la data give us the distance modulus (fi) to each supernova, that is given 
by 

M = m o6s (z i )-M = 51og[— ^] + 25 (3.2) 

Mpc 

where M is their absolute magnitudes. The distance modulus also can be written as 

M = 51og 10 -D L (z) +fi (3.3) 

where Dl = L is the Hubble- free luminosity distance and /io is the zero point offset 
(which is an additional model-independent parameter) defined by 

MO = 5 Iog 10 (— S-) + 25 = 42.38 - 5 log 10 h (3.4) 
Mpc 



In the present paper we used the Union set induing 307 data of Kowalski et. al | |18| ], 
that includes the recent samples from the SNLS []l!J and ESSENCE Surveys |2C( j, older 
datasets, as well as the recently extended dataset of distant supernovae observed with HST 



1 21]. The sample consisting of 414 SNelA was reduced to 307 data after various selection 
cuts were applied in order to create a homogeneous sample. We also used the so called 
"Constitution" set of Hicken et. al |^] induing 397 data, out of which 100 come from the 
new low-z CfA3 sample and the rest from the Union set. Both samples have a redshift 
range of 0.015 < z < 1.55. The main improvement of the Constitution sample is the 
inclusion of a larger number of nearby (z < 0.2) SNela; their inclusion helps to reduce the 



statistical uncertainty [22]. 

The statistic \ 2 1S a useful tool for estimating goodness-of-fit and confidence regions 
on parameters. In our case the XsNia 1S gi ven by 

2 / v _ (fHhe(pi,Zi) - Hobs,i(Zi)) 2 

XSNIaKPi) — 2_^i 2 ^■ 6 - t3 > 

i=l a obs,i 

where pi = (O&0) ^X0> b, B). The x 2 function can be minimized with respect to the hq 
parameter, as it is independent of the data points and the dataset. Expanding the equation 
above with respect to no, we obtain: 
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X 2 (Pi)sNia = A( Pi ) - 2fi B( Pl ) + MoC(Pi) (3-6) 
which has a minimum for fio = B(pi)/Cpi), giving 

2 _ -2 \ / \ B 2 (pi) 

XSNIa,min ~ XSNIa ~ ^{Pi) 777TT V d - ' > 



where 



A(pi) = y (^-^(p^o = o)) 2 (38) 

= £ m-^(p^o = 0) (3 Q) 

Cfo) = ( 3 - 10 ) 
°i 

Now this new XsNia is independent of /iq and can be minimized with respect to the 
parameters of the theoretical model. 

3.2 CMB 

In our investigation we used the method proposed by Komatsu. et al. |2lfl ; in this case 
two distances are important. The distance to recombination can be written as: 



R = V^mOr(zCMB) (3.11) 

where r(zcMB) is given by equation (3.1). The second distance is the angular scale of the 
sound horizon at recombination, r s , that allows to determine the first acoustic peak in the 
following form: 

la = ttI^4 (3.12) 
r s (ZCMB) 



where r, is defined as follows: 



1 f acMB c s (a) 
r s (a C MB) = tt / 2 , . da 3.13 

H Jo a 2 E(a) 



where c s is the sound speed and is given by c s = 1/^/3(1 + 3 ?q b0 ) ■ We are considering 



z cmb = 1090 [23]. Following this reference, we used the prescription for the WMAP5 
distance priors. Thus the \ 2 for the CMB data is: 

XCMB = X CMB C CMB X CMB (3.14) 

where X denote its transpose and the vector X = (R, Ia, 100/i 2 $7bo) and C~ l is the inverse 
covariance matrix, which is given for the WMAP5 observations by: 
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Cij 



/ 0.000367364 0.00181498 -0.000201759 \ 

0.00181498 0.731444 -0.0315874 
^ -0.000201759 -0.0315874 0.00355323 

The best fit turns out to be: 



(3.15) 



XcMB 



R- 1.70 
l a -302.10 

,2 ( 



(3.16) 



\ 100/i^Ofeo - 2.2765 
3.3 BAO 

We used BAO distance measurements obtained at z = 0.20 and z = 0.35 from joint analysis 
of the 2dFGRS and SDSS data [24]. The distance scale used is a combination between the 
radial dilation and the square of the transverse dilation; this is 



Dv(zbao) 



BAO 



dz 



ZBAO 



1/3 



(3.17) 



n H(z) ) H(zbao) 
We also apply the maximum likelihood method and in this case the x 2 1S given by: 



Xbao ~ 

where the inverse covariance matrix is 



r 1 

^BAO 



and 

Xbao r 

where r s is given by equation (3.13). 
3.4 Combining the Datasets 



X BAO C B \ X B AO 



35059 -24031 
-24031 108300 



(3.18) 



(3.19) 



rs(Zdes) 

IV (0.20) 

r S (z d es) 

D v (0.35) 



0.1980 
0.1094 



(3.20) 



We considered that the observational data are independent, so we defined the Xtotal as 



Xtotal 



-2 | 2 | 2 

XsNIa ~r XBAO + Xcmb 



(3.21) 



The best fit values the model can be determined by minimizing the total x 2 - F° r Gaussian 
distributed measurements, the x 2 function is directly related to the maximum likelihood 
estimator. Moreover, if we want to impose a Gaussian prior on one of the parameters being 
measured, pi, centered around pio, with variance dp., we can use the Bayes theorem, and 
write the expression 



L 



cxp 



Xtotal 



exp 



(Pi ~ PiO 

Pi 



,21 



(3.22) 
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Parameter 


Constitution + CMB + BAO 


Union + CMB + BAO 


B 


0.061 ± 0.079 


0.110 ± 0.097 


a 


0.053 ± 0.089 


0.089 ± 0.099 


^xo 


0.834 ± 0.028 


0.854±0.036 


Qo 


-0.60 ± 0.12 


-0.63 ± 0.17 


Zt 


0.62 


0.59 



Table 1: The best fit parameters of the MCG model. The error bars are obtained by marginalized 
likelihood analysis that can be obtained from eq.(3.22). 

In order to constraint the parameters of our interest, we marginalize over the other 
parameters. To account for the uncertainty of the Hubble parameter, we treat it as a free 
parameter and then fix it by using the best fit value of the data. For the reconstruction of 
the q(z) function we used type la Supernovae and BAO data. We followed the standard 
methodology using the Fisher matrix for generating errors. For detais of the method of 
propagation of errors see references |25|. p6| , 

4. Results and Discussion 

The results we obtained for the parameters of the modified Chaplygin equation of state 
P = Bp — Ap~ a are presented in Table 1. We can see that, for both datasets, the exponent 
a is less than 0.1. This shows that the negative pressure component Ap~ a does not differ 
too much from a cosmological constant. The parameter B, related to the dark matter 
component with positive pressure, is also of order 0.1 for both datasets. Theoretically, B 
would be zero for dust-like matter and 1/3 for radiation. The dark energy parameter £lxo 
was found to be approximatelly 0.8, a little bit larger than the 0.7 value of the ACDM 
model. Finally, the kinematical parameter go is of order —0.6 for both datasets, while the 
transition redshift is Zt ~ 0.6. These results for qo and Zt are consistent with constraints 
obtained recently with different methods (H p7|). 

In Figure 2 we show the behavior of the reconstruction of the q(z) function for the 
MGC model. In the cases a, b and c we used the Constitution dataset + BAO data and 
in the cases d, e and / we used the Union dataset + BAO data. In general the results 
using the Constitution sample depend less on the Of,o value. In the case / we can see a 
strong change in the reconstruction of the errors of the q{z) function. In Figure 3 we show 
the one-dimensional likelihood behavior in the case Constitution + BAO + WMAP5 for 
the MCG model. In the upper panels of Figure 4 we show observational constraints using 
only type la Supernovae data. These results change with the inclusion of the BAO data, 
but if we include the WMAP5 data we obtain results which strongly reduce the space of 
parameters. These results are displayed in Figure 5, in which we clearly see that for both 
samples the a = value is included within the confidence region of the space parameters. 
The ACDM model is the particular case a = and B = 0. Therefore, the points with 
a = in Figure 5 do not represent the ACDM model, because we used the best fit for the 
B value (see Table 1); the finite value of B is the distinctive feature of the MCG model. 
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Figure 2: In panels a,b,c we present the best fit and reconstruct the erros for the Constitution 
datasct + BAO, using B = 0.061, and in panels d,e,f for the Union datasct + BAO, using B = 0.110. 
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Figure 4: In the top left we show the confidence regions for the Constitution dataset and in the 
top right for the Union dataset. In the bottom left confidence regions for the case of Constitution 
+ BAO and in the bottom right for the case Union + BAO. In all cases we used B given by Table 
1 and a prior of il^ = 0.042. 
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Figure 5: In the top left we show observational constraints of SNela Constitution + BAO + CMB. 
In the top right, SNela Union + BAO + CMB. Each bottom panel shows the intersection of the 
curves above. 
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